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1.  Introduction 

The  major  aim  of  our  research  investigations  has  been  to  develop  new  algorithms  and  the 
corresponding  computer  codes  for  exact  numerical  solutions  of  the  time-dependent  Boltzmann 
equation  for  the  transport  of  electrons  (and  other  charged  particles)  in  gases  subjected  to  external 
electric  fields  of  arbitrary  strengths.  The  electron  velocity  distribution  function  obtained  as  a 
solution  of  the  Boltzmann  equation  is  used  to  calculate  a  set  of  experimentally  measurable  swarm 
parameters  for  the  transport  of  electrons  in  actual  pure  gases  or  in  gas  mixtures  subjected  to  static 
or  time-dependent  external  electric  fields.  A  direct  comparison  of  the  calculated  swarm  parameters 
with  the  corresponding  experimental  values  provides  a  useful  check  on  the  validity  and  accuracy  of 
the  results  obtained  using  the  present  codes.  It  is  especially  useful  to  have  a  general  computer 
program  for  calculating  the  theoretical  values  of  the  electron  swarm  parameters  for  gases  for  which 
there  are,  as  yet,  no  experimental  values  or  for  which  direct  laboratory  experiments  are  simply  not 
feasible. 

Using  our  novel  algorithm  we  have  obtained  the  numerical  values  of  swarm  parameters  for 
the  cases  in  which  the  external  electric  field  is  either  constant  or  varies  with  time  in  a  sinusoidal 
fashion,  both  for  a  pure  gas  as  well  as  for  a  gas  mixture.  The  external  rf  electric  field  is  of  the 
form  E  =  Eq  sin  cot  or  E  =  E^^  cos  cot,  of  arbitrary  strength  and  arbitrary  frequency  co.  The 

present  investigations  have  allowed  us  to  analyze  the  relations  between  the  frequency  (O  of  the  field 
and  the  frequency  v  of  collisions  of  the  electrons  with  background  gas  particles,  with  number 
density  N.  These  relationships  can  be  qualitatively  understood  by  using  a  simple  model  in  which 
one  assumes  a  constant  collision  frequency  for  the  projectile-gas  collisions.  In  this  model  the  time 
dependence  of  various  swarm  parameters  using  the  complete  Boltzmann  equation  can  be 
analytically  solved.  The  whole  analysis  is  carried  out  for  an  arbitrary  value  of  the  ratio  of  the 
projectile  mass  to  the  ambient  gas  particle  mass  so  that  the  derived  expressions  are  valid  even  for 
the  transport  of  heavy  ion  swarms  in  gases  subjected  to  external  rf  fields  of  arbitrary  strength. 
Even  though  these  analytical  results  are  interesting  and  useful  in  their  own  right,  they  have  been 
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used  merely  as  a  check  on  the  results  from  the  actual  computer  code  which  did  not,  in  general, 
assume  a  constant  collision  frequency;  in  the  code  the  energy-dependent  collision  frequency  is 
determined  by  the  actual  collision  cross  sections  used. 

In  the  course  of  present  work,  we  came  across  a  few  additional  problems  that  needed  to  be 
investigated.  For  example,  in  order  to  test  the  versatility  of  our  single-target-gas  computer  code  we 
tried  to  obtain  the  temporal  history  of  the  behavior  of  charged  particles  other  than  electrons  in  a 
gas,  and  for  this  purpose  we  chose  a  positron  swarm.  This  investigation  allowed  us  to  compare 
the  behavior  of  electron  swarms  with  that  of  positron  swarms  under  similar  conditions.  As  a  part 
of  additional  investigations,  we  calculated  the  cross  sections  for  the  processes  which  are  important 

for  a  realistic  modeling  of  a  hydrogen  plasma.  In  this  connection  we  have  investigated  the 
dependence  of  the  cross  sections  for  pure  dissociation  of  H2  upon  the  initial  vibrational  state  of  H2. 

Finally,  we  have  also  investigated  the  penetration  probability  of  a  quantum  harmonic  oscillator  into 
the  classically  forbidden  region  in  the  limit  of  large  quantum  numbers  using  analytical  methods.  An 
expression  is  derived  for  the  asymptotic  form  of  the  penetration  probability,  and  is  shown  to  be  in 
good  agreement  with  an  earlier  result  obtained  by  numerical  methods.  The  asymptotic  form  of  the 
probability  density  in  the  region  between  the  classical  turning  points  is  also  presented  and  found  to 
have  a  simple  physical  interpretation. 

1.1  Novel  method  of  solving  the  Boltzmann  equation 

We  provide  here  a  very  bnef  summary  of  our  algorithm  for  solving  the  Boltzmann  equation 
exactly.  The  traditional  Boltzmann  equation  used  for  calculating  the  parameters  which  characterize 
the  transport  of  electrons  injected  in  a  neutral  gas  which  is  subjected  to  an  external  electric  field  E 
(constant  or  time-dependent)  is 

+  a  »  Vy  f(v,  t)  =  R(v,t)  ,  (1) 

The  solution  f(v,  t)  provides  the  (spatially  homogeneous)  velocity,  or  equivalently  energy, 
distribution  function  of  the  electrons.  Here  a  =  -eE/m  represents  the  acceleration  of  the  electrons 
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due  to  external  electric  field  E  and  the  collision  term  R(v,  t)  represents  the  rate  of  change  in  the 
electron  velocity  distribution  function  due  to  all  possible  collision  processes  among  the  electrons 
and  the  ambient  gas  particles.  The  collision  term  R(v,  t)  itself  involves  integrals  over  the 
distribution  function  f(v,  t)  so  that  in  its  complete  form,  the  Boltzmann  equation  is  a  time- 
dependent  integrodifferential  equation  whose  numerical  solution,  in  the  past,  has  been  obtained 
only  after  making  some  simplifying  assumptions. 

In  our  algorithm  no  expansion  of  any  kind  of  the  distribution  function  is  made  and  no 
derivatives  are  evaluated  numerically.  The  proposed  procedure  for  obtaining  the  time-dependent 
velocity  distribution  function  is  numerically  stable  and  computationally  inexpensive.  In  order  to 
derive  the  basic  equation  of  our  algorithm  we  multiply  Eq.  (1)  by  a  small  time  interval  At  and  add 

the  distribution  function  f(v,  t)  to  it  to  obtain 

f(v,  t)  +  ^At  +  Av  •  Vy  j  f(v,  t)  =  f(v,  t)  +  R(v,  t)  At  ,  (2) 

where  Av  =  a  At  is  a  small  velocity  increment.  Eq.  (2)  can  be  rewritten  as 

f(v  +  Av,  t  +  At)  =  f(v,  t)  +  R(v,  t)  At  .  (3) 

Thus,  the  difference  equation  (3)  is  entirely  equivalent  to  the  Boltzmann  equation  in  its  physical 
content.  We  have  realized  [1]  that  the  above  difference  equation  is  much  simpler  to  solve  for  the 
time-dependent  velocity  distribution  function  than  the  standard  integrodifferential  Boltzmann 
equation.  Furthermore,  being  a  difference  equation,  Eq.  (3)  is  more  suited  for  numerical 
computations  than  Eq.  (1).  A  few  important  advantages  of  our  novel  algorithm  over  the  traditional 
techniques  are: 

(i)  Since  our  algorithm  does  not  involve  any  numerical  computation  of  partial  derivatives,  there 
is  no  need  to  invoke  Courant-Levy  condition  for  numerical  stability.  Thus,  the  time-step  At 
and  the  velocity- step  Av  can  be  varied  almost  independently. 
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(ii)  Lack  of  numerical  evaluation  of  any  derivatives  also  makes  our  algorithm  computationally 
inexpensive.  Almost  all  of  our  numerical  results  are  obtained  on  microcomputers  or 
workstations. 

(iii)  Since  no  truncation  of  expansion  of  any  kind  of  the  distribution  function  is  made  in  our 
work,  the  present  algorithm  provides  swarm  parameters  for  any  value  (that  is,  for  both  high 
and  low  values)  of  E/N. 

1.2  The  collision  term 

An  explicit  expression  for  the  collision  term  R(v,  t)  has  been  provided  by  Holstein  [2]  in 
the  following  form: 

R(v,t)  =  5^R^(v,t)  -  R:;.(v,t)  ,  (4) 

P 

where 

N  f  ^ 

^  J  '"p  0  J  ''p  '  Sp(Vp,  V))  (5a) 

represents  the  rate  at  which  the  projectile  particles  with  initial  speed  Vp  are  scattered,  due  to  the  pth 
(p  =  elastic,  excitation,  ionization  etc.)  scattering  process,  into  a  velocity-space-volume  element 

ri 

d  V  located  at  v,  and 

R^(v,  t)  =  N  V  f(v,  t)  CT.j.(v)  (5b) 

represents  the  rate  at  which  the  projectile  particles  are  scattered,  due  to  aU  possible  collision 
processes,  out  of  the  velocity-space-volume  element  d\  located  at  v.  Naturally  <^p(Vp,  v)  is  the 

differential  cross  section  for  the  pth  collision  process,  and  <Tj(v)  is  the  integrated  total  scattering 
cross  section.  The  energy  conserving  delta-function  relates  the  initial  speed  Vp  to  the  final  speed  v 
for  the  pth  collision  process  via  the  function  gp(Vp,  xja)  which  is  defined  by  the  relationship 

V  =  gp(Vp,  tk). 
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V  is  the  scattering  angle.  Figure  1  shows  the  relationship  among  various  angles  as  well  as 
directions  of  relevant  vector  quantities.  For  example,  the  direction  of  Vp  can  be  specified  with 

respect  to  either  v  or  z.  These  two  specifications  imply 

d^Vp  =  Vp  dVp  sin  0p  d0p  d(j)p  =  Vp  dVp  sin  \it  d\|/  da  .  (6a) 

Also,  cos  \|/  =  cos  0  cos  0p  +  sin  0  sin  0p  cos  (l)p  .  (6b) 


Fig.  1:  Geometry  used  in  the  discussion  of  the  collision  term. 


For  elastic  collisions  (p  =  e)  of  projectiles  of  mass  m  incident  on  gas  particles  (of  mass  M) 
at  rest,  we  obtain  from  simple  kinematics  (simply  using  conservation  of  energy  and  momentum), 

(7) 
(8a) 


V  =  Vgq(\|/)  s  g^(v^,\|/) 


Here 


e'  e’ 

_  ( 1  -  sin^\|t)  +  q  cos  V 
1  +  p. 


q(v)  = 

and  p  =  m/M.  For  small  p  (for  example,  for  electron  or  positron  swarms  in  atomic  and  molecular 
gases)  we  can  approximate 


q(v)  =  1  -  P  (1  -  cos  v) 


(8b) 
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Now,  no  matter  what  simplifying  assumptions  one  makes  for  obtaining  a  solution  of  Eq. 
(1),  ultimately  one  will  compute,  for  direct  comparison  with  the  corresponding  experimental 
values,  various  transport  or  swarm  parameters  (such  as  the  average  energy  of  the  projectile 
particles  in  the  charged  beam,  the  drift  and  diffusion  rates  of  the  beam  particles  in  the  ambient  gas, 
the  rates  of  excitation  and  ionization  of  the  ambient  gas  particles  by  the  impact  of  the  charged  beam 
particles  etc.),  which  are  expressed  as  integrals  involving  the  distribution  function  f(v,  t).  A 
typical  swarm  parameter  will  be  the  expectation  value  of  a  function  g(v)  =  g(v,  0,  (j)),  which  itself 


does  not  depend  on  time  explicitly: 


G(t)  -  <g(v))  = 


(9) 


f  2 

where  n(t)  =  2k  J  dv  /  f(v,  0,  t)  sin  0  d0  is  the  number  density  of  projectiles  and,  in  the 
0  0 


case  of  elastic  collisions  only,  n(t)  is  independent  of  time  and  can  be  normalized  to  1  for 
convenience.  Note  that  f(v,  0,  t)  is  azimuthally  symmetric  since  E  will  be  taken  to  be  along  the  z 


axis. 


ision 


Because  of  the  potential  applications  of  a  hydrocarbon-rare-gas  mixture,  we  chose  to 

investigate  the  transport  of  electrons  in  a  mixture  of  gaseous  methane  and  argon  with  the  ratio  of 
the  number  densities  of  CH^  and  Ar  always  maintained  at  0.05.  In  the  collision  term  of  the 


Boltzmann  equation  we  included  the  cross  sections  for  elastic,  vibrational  excitation  and  ionization 
processes  for  methane,  and  for  the  elastic,  excitation  and  ionization  processes  for  argon.  Using  the 
numerical  electron  velocity  distribution  function  we  calculated,  as  a  function  of  time,  various 
swarm  parameters  such  as  the  transverse  diffusion  coefficient,  drift  velocity  and  average  energy  of 
the  electron  swarm,  rates  of  ionization  and  excitation  of  the  background  gas  mixture  by  electron 
impact,  frequencies  of  collisions  of  electrons  with  ambient  gas  targets,  etc.  For  brevity  of 
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presentation,  we  will  only  show  the  time-dependence  of  a  single  swarm  parameter,  namely,  the 
average  energy  of  electrons  in  the  swarm.  In  all  presentations,  the  initial  (at  t  =  0)  velocity 

distribution  of  the  electrons  in  the  swarm  is  assumed  to  be  a  Maxwellian  distribution  with  an 
average  energy  of  10  eV.  Also,  at  t  =  0,  various  collision  frequencies,  defined  as  Vp  =  <NapV> 

(p  =  elastic,  excitation,  ionization,  etc.),  are  calculated  for  both  methane  and  argon  using  the 
Maxwellian  distribution  function  for  electrons.  The  largest  and  the  most  relevant  of  these 
frequencies,  in  the  present  mixture,  is  the  frequency  of  elastic  collisions  of  electrons  with  argon. 
Note  that  the  electron  velocity  distribution  function,  and  consequently  the  collision  frequencies, 
varies  with  time;  however,  at  all  times  the  frequency  of  elastic  collisions  of  electrons  with  argon 
remains  as  the  dominant  collision  frequency.  In  the  present  investigations,  we  have  varied  co  while 
keeping  v  fixed.  The  results  of  these  investigations,  shown  in  fig.  2,  provide  insight  in  the  time- 
dependent  behavior  of  electron  swarms  in  gases  subjected  to  external  rf  electric  fields. 

Fig.  2,  consisting  of  five  frames  labeled  (a)  to  (e),  shows  the  temporal  evolution  of  the 

average  electron  energy  in  the  methane-argon  mixture  with  the  value  of  the  field  frequency  O) 
increasing  from  fig.  2(a)  to  fig.  2(e).  The  set  of  parameters  belonging  to  this  figure  are:  =  744 

V/cm,  the  total  gas  density  N  =  3.72  x  10^^  cm*^,  =  2  kTd,  field  frequency  O)  =  1.0  x  10^ 

Hz,  v^^  =  5.80  X  10^  sec*^  N(CH4)  =  1.77  x  10^^  cm'^,  N(Ar)  =  3.54  x  10^^  cm'^  and 
N(CH4):N(Ar)  =  1 :20.  On  advancing  from  fig.  2(a)  to  fig.  2(e),  the  electric  field  frequency  to  is 
increased  while  keeping  E^,  N  and  the  ratio  N(CH4):N(Ar)  fixed.  In  figs.  2(a)  and  2(b)  the 

collision  frequency  v  is  larger  than  the  field  frequency  co;  in  fig.  2(c)  the  collision  frequency  v  is 
comparable  to  the  field  frequency  to  and  in  figs.  2(d)  and  2(e)  the  collision  frequency  v  is  less  than 
the  field  frequency  to.  Results  in  fig.  2  allow  us  to  make  the  following  conclusions: 

•  We  observe,  from  figs.  2(a)  and  2(b),  that  when  collision  frequency  v  is  larger  than  the  field 
frequency  co,  the  average  of  the  oscillating  equilibrium  values  of  the  electron  energy  is  more 
than  the  initial  electron  energy  (which,  in  the  present  case,  is  10  eV  at  t  =  0)  leading  to  the 


Average  Electron  Energy  (eV) 
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Time  (sec) 

Fig.  2:  Temporal  evolution  of  the  average  energy  of  an  electron  swarm  in  a  methane-argon  mixture 
subjected  to  an  external  if  electric  field  for  different  values  of  the  field  frequency,  O). 
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microwave  heating  of  the  electron  swarm.  For  v  less  than  ©  (see  figs.  2(d)  and  2(e))  the 
electron  swarm  loses  energy  to  the  background  gas  particles. 

•  As  the  frequency  ©  of  the  field  is  increased,  the  amplitude  of  the  oscillations  in  the  average 
electron  energy  is  reduced  (compare  figs.  2(a)  through  2(e)).  Furthermore,  for  very  high 
values  of  ©,  the  average  electron  energy  behaves  (see  fig.  2(e))  in  the  same  manner  as  in  the 
absence  of  external  field.  Both  of  these  facts  can  be  physically  understood  by  realizing  the 
origin  of  oscillations  in  the  time-dependence  of  the  average  electron  energy.  Change  in  the 
average  energy  of  the  electrons  is  caused  by  collisions  as  well  as  by  the  external  electric  field; 
the  external  rf  field,  during  its  one  full  cycle,  both  accelerates  as  well  as  decelerates  the 
electrons  so  that  the  average  energy  of  electrons  exhibits  the  oscillatory  behavior.  As  the  field 
frequency  ©  increases  (and,  therefore,  the  period  decreases)  the  time  available  for  the  electrons 
to  be  accelerated  or  decelerated  by  one  cycle  of  the  field  is  reduced  which  leads  to  a  diminished 
amplitude  of  the  oscillation  in  the  time-dependence  of  the  average  electron  energy.  At 
extremely  high  frequency  of  the  field  (or  very  small  period  of  oscillation  of  the  field),  a  typical 
electron  makes  almost  no  collisions  during  one  period  and,  furthermore,  both  acceleration  and 
deceleration  of  the  electron  by  the  field  cause  no  net  change  in  the  velocity  of  the  electron  -  any 
change  in  the  velocity  of  the  electron  occurs  due  to  collisions  only.  Thus,  in  fig.  2(e)  where  © 
is  much  larger  than  v  the  average  electron  energy  behaves  as  if  there  is  no  external  field.  A 
physical  understanding  of  these  results  can  be  obtained  using  a  simple  model  (see  sec.  1.4) 

The  investigations  of  electron  transport  in  a  gas  mixture  have  demonstrated  that  the  electron 
swarm  parameters  depend  sensitively  on  the  relative  values  of  the  electron  collision  frequency  v 
and  the  external  field  frequency  ©.  Such  investigations  are  important  in  developing  a  tailor-made 
mixture  of  gases  with  a  pre-chosen  behavior  for  electron  transport  in  them. 

1.4  Constant  Collision  Frequency  Model. 

While  developing  our  computer  code  for  the  case  of  a  sinusoidally  varying  time-dependent 
electric  field,  E  =  z  cos  (©t),  we  intended  to  consider  simple  cross  sections  in  the  collision 
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term  such  that  the  complete  time-dependence  of  the  corresponding  swarm  parameters,  G(t),  can  be 
exactly  obtained  analytically.  These  analytical  expressions  for  the  time-dependent  swarm 
parameters  will,  then,  provide  useful  checks  on  the  correctness  of  the  numerical  results  obtained 
from  our  general  computer  code  with  time-dependent  electric  field.  With  this  thought  in  mind  we 
started  with  the  case  where  one  only  has  isotropic  elastic  scattering  of  electrons,  positrons  or  heavy 

ions  by  atomic  and  molecular  targets  with  a  constant  frequency  v.  Then 
V  =  N  V  a.j.(v)  =  4'7t  N  v  ct^(v,  xj/) 


Jt 


2  K 

i' 


Vg  dVg  j  sin  V  dv  f  da  f{Vg,  t)  5(v  -  v^  q(\jr)) 


and 


R'Cv,  t)  =  v  f(v,  t) 


(10) 

(lla) 

(llb) 


Now  we  can  write  the  differential  equation  satisfied  by  the  swarm  parameter  G(t)  [with  n(t) 


normalized  to  1]  simply  by  differentiating  Eq.  (9), 

2  7t 

dG(t) 

—  J  V  uv 

0 


dt 


0 

af 


TZ 

J  v^  dv  J  sin  0  d9  J  d(t)  g(v,  0,  <[))  • 

)  0  ^ 


(12) 


Substituting  for  ^  from  the  Boltzmann  equation  (1),  taking  a  =  qE/m  to  be  along  the  z-axis  and 


using  integration  by  parts,  we  obtain 


dt  ■  ^  ^ 


ag(v,  0,  0)  sin  0 
a0  V 


-  cos 


Q  ag(v,  0, 0) 

5v 


>  +  V 


dGjj^(t) 

dt 


-  (g(v,  0,  0)>  ,  (13) 


where,  as  before,  (...)  represents  the  expectation  value  and 

_  2n 

dGiR(t) 


dt 


TC 

Jv^dv  J  sin  0  d0  /  d0  g(v,  0,  0)  R^(v,  t) 


(14) 


Next,  substituting  R  (v,  t)  from  Eq.  (11a)  into  Eq.  (14),  using  (6a)  and  using  the  delta  function  to 


carry  out  the  v-integration,  we  obtain 

oo 


2  K 
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where  the  function  H  is  given  by 

71 


-  71  27t 

H(Vg,  e^,  <t)g)  =  4^  J  sin  e  d9  j  d(t)  g(v^  q(\|7),  0, 0) 


(16) 

0  0-' 

and  w  is  related  to  0,  0  and  6  via  Eq.  (6b).  Thus,  we  obtain  our  principal  result,  namely, 

the  first-order  ordinary  differential  equation  satisfied  by  a  swarm  parameter  G(t) : 
dG(t)  _  d(g(v,  0,  0)) 
dt  dt 


=  -a< 


cos 


)  +  V  (H(v,  0,  (j))  -  g(v,  0,  ([))>  ,  (17) 


where  H(v,  0,  ({))  is  given  by  equation  (16).  As  we  will  see  now,  there  are  several  important  cases 

in  which  the  differential  equation  (17)  can  be  solved  exactly.  In  all  these  cases  we  will  take  the 
external  electric  field  to  be  sinusoidally  time-dependent,  that  is,  E  =  z  E^  cos  (wt),  or 
equivalently,  a  =  a^  cos  (tot).  We  consider  two  specific  examples: 

fa )  Drift  velocity  of  electrons. 

In  this  case  g(v,  0,  ([))  =  v  cos  0,  and  <g(v,  0,  ([)))  =  v^(t),  the  drift  velocity.  Using  Eqs. 
(16),  (6b)  and  (8)  together,  we  get,  in  this  case, 

.  3  (1^.  \  (6,)- 

The  differential  equation  (17),  now,  becomes 

+  a^  cos  (cot)  -  vv^(t)  ,  (18) 

where  v  =  (1  -  X)  v  and  ^  ^  differential  equation  can  be  solved  exactly;  the 

solution  is 

Vd(t)  =  Vd(0)  exp {-  V  t}  -  A  [  V  {exp(-  v  t)  -  cos  (cot) }  -  to  sin  (tot)],  (19a) 

with  A  =  .  (19b) 

V  -I-  to 

As  a  special  case,  for  a  dc  electric  field  (to  =  0),  one  has 

Vd(t)  =  Vd(0)  exp{-  V  t}  +  ^  [1  -  exp  (-  v  t)]  .  (20) 


In  the  absence  of  collisions  (v  =  0)  Eq.  (20)  reduces  to 

Vd(t)=Vd(0)  +  a^t  , 


(21) 
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which  is  consistent  with  the  simple  equation  of  kinematics. 

fb)  Average  energy  of  electrons. 

For  the  case  of  average  energy  of  electrons,  g(v,  0,  q)  =  2^'^ 


For  this  case,  we  obtain  from  Eq.  (16)  [using  (6b)  and  (8)], 

n 


.  1  2:1  ^ 

H(Vg,  6^,  (t)g)  =  jmvg  ^  I  d(l)  J  qV)sin0d0  = 
and  <H(v,  0,  <|))>  = 


.3(l+|i)^J 


’^^le(t). 

.3  (1 + ]xn 


The  differential  equation  for  e(t)  now  becomes 


d£(t) 

dt 


=  +  m  a„  cos  (cot)  v  .(t)  +  v  \ — ^  ~  ^  ^  £(0 

o  ^  ^  l3(l+p)2  J 


or 


^4^  =  +  m  a^  cos  (cot)  v  .(t)  -  6  X  v  e(t) 


dt 

This  differential  equation  for  e(t)  can  be  solved  exactly;  the  solution  is 

t 


e(t)  =  exp  {-  6  ^  V  t} 


e(0)  +  ni  Jdt'  cos  (cot')  exp  {6  X  v  t’)  v^(t') 


(22a) 


(22b) 


(23) 


After  substituting  v^(t)  from  Eq.  (19)  into  Eq.  (23)  one  can  easily  carry  out  the  integrals  and  the 
resulting  expression  for  e(t)  is 

e(t)  =  exp{-6Xvt}  [e(0)  +  m 3^1(0],  (24a) 

where 


V  .(0)  -  A  V  A 

I(t)  =  ■-ii.i-z — 2  'W  +  7 

(1 -6A)^v^  +  to^  ^ 


■-  3X^  +  co2' 


9^2  7^2.  _L_  ,.,2 


'  V""  +  0" 


+  6Xvt}  -  1] 


+ 


A 

4 


~exp  {  6  X  v  t)' 
.  9X^  ^  . 


[cos  (2wt){3X^  -  ©2)  +'  sin  (2©t)  {(1  +  3X)  v  ©}] 


with 

I(t)  =  exp  {-  (1  -  6X)  V  t}  [©  sin  (©t)  -  (1  -  6X)  v  cos  (©t)]  +  (1  -  6X)  v  .  (24b) 

Again,  for  the  case  of  a  dc  external  field  (©  =  0),  Eq.  (24)  simplifies  to 
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1 

Time  (seconds) 

Fig.  3:  Time-dependence  of  the  electron  drift  velocity  in  a  model  gas,  with  constant  collision 
frequency  v,  subjected  to  an  external  dc  electric  field  (ca  =  0). 


Collision  Frequency,  v  =  20  GHz 

-  Numerical  Results 

-  Analytical  Results 


Time  (seconds) 


Fig.  4:  Time-dependence  of  the  electron  average  energy  in  a  model  gas,  with  constant  collision 
frequency  v,  subjected  to  an  external  dc  electric  field  (©  =  0). 
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xlO^ 


Fig.  5:  Time-dependence  of  the  electron  drift  velocity  in  a  model  gas,  with  constant  collision 
frequency  v,  subjected  to  an  external  ac  electric  field  (frequency  to). 


Fig.  6:  Time-dependence  of  the  electron  average  energy  in  a  model  gas,  with  constant  collision 
frequency  v,  subjected  to  an  external  ac  electric  field  (frequency  to). 
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e(t)  =  exp  {-  6  ^  V  t}  e(0)  +  ^  \  {-  6  X  v  t}  -  exp  {-  v  t}] 

(1  -  oA)  V 


m  a, 


:^[l  -  exp  {-6A,vt}] 


6Xv 

Furthermore,  in  the  absence  of  collisions  (v  =  0),  Eq.  (25)  becomes 

majt^ 

e(t)  =  £(0)  +  m  v^(0)  t  +  — ^ 


(25) 


(26) 


which  is  consistent  with  Eq.  (21),  the  simple  equation  of  kinematics.  Considerable  simplification 
is  also  achieved  in  Eq.  (24)  for  the  case  where  |a  =  0;  this  is  a  quite  reasonable  assumption  when 
electron  or  positron  swarms  are  transported  in  atomic  and  molecular  gases  since  then  the  ratio  m/M 
is  indeed  quite  small.  In  this  case  one  obtains,  on  setting  A.  =  0, 


e(t)  =  e(0)  +  m  A  [v^(0)  -  A  v]  I(t)  +  m  A  a^ 


1  -  cos  (2cot)  V 
4 


sin  (2(ot)  y_tl 
4  CO  “^2]’ 


where  I(t)  =  exp  (-  v  t)  [co  sin  (cot)  -  v  cos  (cot)]  +  v  and  A  =  ' 

Equations  (19),  (25)  and  (27)  have  provided  very  useful  checks  on  the  computer  code  which 
incorporates  a  sinusoidally  time-dependent  external  electric  field.  For  example.  Figs.  3  and  4 
show  a  comparison  of  the  time-dependence  of  the  electron  drift  velocity  and  electron  average 
energy,  in  a  model  gas  with  constant  collision  frequency  and  subjected  to  an  external  dc  electric 
field,  obtained  from  the  numerical  code  and  from  the  present  analytical  work  [see  Eqs.  (20)  and 
(25)].  Figs.  5  and  6  show  a  similar  comparison  for  the  case  in  which  the  external  electric  field  is 
sinusoidally  time-dependent  with  a  frequency  ca  [see  Eqs.  (19)  and  (24)].  These  comparisons  have 
not  only  provided  useful  checks  on  the  numerical  results  for  electrons  in  real  gases  but  have  also 
enabled  us  to  decide  the  actual  time-step  At  and  velocity  step  Av  which  should  be  used  for  accurate 
numerical  results,  for  both  ac  and  dc  electric  fields.  Eq.  (27)  can  be  used  to  understand,  at  least 
qualitatively,  the  conclusions  reached  from  the  results  of  Fig.  2. 
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1.5  Effect  of  the  temperature  of  the  background  gas  on  f(v.  tl. 

In  the  case  where  the  background  gas  particles  are  at  rest,  even  for  elastic  collisions  of 
electrons  (of  mass  m)  with  gas  particles  (of  mass  M)  there  is  a  transfer  of  energy,  of  the  order  of  |i 
=  m/M,  from  the  electron  to  the  recoiling  gas  particle  so  that,  in  the  absence  of  any  external  electric 
field,  the  numerical  value  of  the  average  energy  of  electrons  continues  to  reduce  with  time, 
approaching  the  unrealistic  value  of  zero.  A  finite  equilibrium  value  of  the  average  electron  energy 
is  reached,  numerically,  only  if  one  assumes  |i  =  0.  In  order  to  make  the  situation  in  our 
calculations  more  realistic,  we  decided  to  investigate  the  effect  of  temperature  (and,  therefore, 
motion)  of  the  background  gas  on  the  electron  velocity  distribution  function.  We  assume  that  the 
number  density  of  the  gas  particles  is  much  larger  than  that  of  electrons  so  that  electrons  interact 
with  gas  particles  more  frequently  than  among  themselves.  Also,  because  of  their  large 
concentration,  the  gas  particles  interact  with  each  other  much  more  often  than  with  electrons. 
Therefore,  the  background  gas  can  be  assumed  to  have  a  Maxwellian  velocity  distribution  function, 
with  a  temperature  parameter  that  varies  slightly  with  time.  The  collision  term,  in  general,  has  the 
form 

R(v,  t)  =  J{  f(v'.  t)  F(V’,  t)  -  f(v,  t)  F(V,  t) }  Iv  -  VI  dQ  dV  ,  (28) 


where  v'  and  V  are  the  velocities  of  the  electron  and  the  gas  particle  before  the  collision  and  v 
and  V  are  the  corresponding  velocities  after  the  collision.  F(V,  t)  is  the  Maxwellian  velocity 


distribution  function  of  the  gas,  Iv  -  VI  =  u  is  the  relative  velocity  of  the  colliding  partners  and  ^ 
is  the  differential  cross  section  for  elastic  scattering.  In  cylindrical  coordinates  the  collision 


integral,  for  small  but  finite  p,  can  be  written  as 

71  2jt  7t  271 


^11  II  . 


N  C!^  f(v  v^  ,  t)  V 


(29) 
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where,  as  shown  in  fig.  7,  and  a  are  the  polar  and  azimuthal  angles  of  V,  (0,  <|))  are  the 
scattering  angles,  is  the  integrated  elastic  cross-section  and  N  is  the  gas  number  density.  The 

five-dimensional  collision  integral  is  evaluated,  in  the  present  work,  numerically  using  the  Monte 
Carlo  technique.  To  do  it  we  need  to  express  v  ',  v  '  and  V  in  terms  of  v  ,  v  ,  V  and  all  the 

p  z  P  z 

relevant  angles  (\|/,  a,  0  and  (j))  using  various  conservation  laws. 

From  the  energy  and  momentum  conservation  laws  we  know  that,  in  elastic  collisions,  the 
vector  of  center  of  mass  velocity  and  the  magnitude  of  relative  velocity  of  the  two  colliding  paitners 
do  not  change.  The  center  of  mass  velocity  and  relative  velocity,  given  by 


C  = 


mv-HMV  mv'-i-MV' 


m-i-M 


m  +  M 


,  u  =  V  -  V  ,  u'  =  V  -  v',  [with  u  =  u'] 


have  components 

S' 


V 

p 


V  2l‘''pVpC0sa  ,  q  .  T4Tr<t‘’'z+S)’ 


“p  =  V  ''p^  +  V  ■  ^''p  U^  =  v^-V^,and  u  =  V  UpSu/  . 

These  relations  provide  the  components  of  the  center  of  mass  velocity  and  the  magnitude  of  relative 

velocity  in  terms  of  the  velocities  of  particles  after  the  collision.  Now  we  need  to  obtain  them  in 

terms  of  the  velocities  before  the  collisions.  Using 
V  =  C  -  ^  u’  , 


l  +  ii 


the  components  of  v'  and  of  u',  in  terms  of  scattering  angles,  are 

V'-n/ 


c}+ 


p  /I  J.  ■1^2  p  1  IJ.  “p  '^p 


(i  +  iir 


u  '  C„  cos  d) 

n  n  ' 


''z’  =  S-TT]I“z’. 


with  Up'  =  u  sin  0  ,  =  u  cos  0. 

Various  expressions,  given  above,  combine  together  to  provide  the  velocity  of  electron  before  the 
collision  in  terms  of  the  required  variables.  Finally,  we  get  V  using  the  energy  conservation  law: 
V’  =  V  lt(v^- v’^)-f-V^  . 

Now  we  can  carry  out  the  five-dimensional  integral  using  the  standard  Monte  Carlo  technique,  that 
is,  by  generating  random  values  of  V,  0,  \jr,  (J),  and  a  and  by  calculating  the  corresponding 
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integrand.  Using  this  collision  integral,  which  incorporates  the  effect  of  temperature  of  the 

background  gas,  in  the  Boltzmann  equation  we  can  evaluate  the  electron  velocity  distribution 

function.  For  this  purpose,  we  investigated  distribution  of  electrons  in  gaseous  argon.  The  initial 

2  2 

distribution  of  electrons  is  assumed  to  be  a  Maxwellian  function,  ^th)’ 

related  to  the  electron  temperature  T^  (which  is  taken  in  this  model  investigation  as  10^  K)  and  the 
initial  gas  temperature  is  taken  to  be  one-half  of  the  electron  temperature.  For  the  mass  ratio  \i  = 
10'^,  which  is  of  correct  order  of  magnitude  for  electrons  in  argon,  fig.  8a  compares  the 

numerically  calculated  equilibrium  electron  velocity  distribution  function  (solid  line)  with  the  initial 
Maxwellian  distribution  function  (dashed  line)  for  Vp  =  0  and  for  no  external  electric  field. 


-2  0  2 
Vz/Vth 


Fig.  8:  Electron  velocity  distribution  function  for  two  different  values  of  the  mass  ratio,  |i. 
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Obviously,  the  effect  of  temperature  of  the  ambient  gas  is  hardly  noticeable  in  this  figure. 
However,  if  the  magnitude  of  the  mass  ratio  |a  is  artificially  inflated  to  an  unrealistic  value  of  10’  , 
the  electron  velocity  distribution  function,  at  equilibrium,  differs  perceptibly,  as  seen  in  fig.  8b, 
from  the  initial  Maxwellian  function.  Based  on  the  results  of  fig.  8  it  appears  that  the  effect  of 
temperature  of  the  background  gas  on  the  electron  velocity  distribution  function,  and  consequently 

on  the  swarm  parameters,  may  not  be  important  as  long  as  there  is  no  external  electric  field. 

2.  Comparative  study  of  electron  and  positron  transport  in  gases. 

In  order  to  compare  the  time-dependent  behavior  of  a  swarm  of  electrons  in  various  atomic 
gases  with  the  behavior  of  swarms  of  other  charged  projectiles,  we  have  investigated  the  transport 
properties  of  electrons  and  positrons  in  the  same  ambient  gas.  Positrons,  being  antiparticles  of 
electrons,  are  quite  similar  to  electrons  in  all  respects  except  that,  unlike  electrons,  the  polarity  of 
their  charge  is  positive.  Investigations  of  the  comparative  behavior  are  accomplished  by  solving 
the  Boltzmann  transport  equation  exactly  using  our  novel  algorithm  which  has  been  described 
above.  For  both  electrons  and  positrons,  we  have  to  include,  in  the  collision  term  of  the 
Boltzmann  equation,  the  cross  sections  for  the  elastic,  excitation  and  ionization  processes.  In 
addition,  in  the  case  of  positrons,  the  effect  of  positronium  formation  and  of  positron  annihilation 
had  to  be  included  to  achieve  the  correct  equilibrium  values  of  various  swarm  parameters  (such  as 
the  average  energy  and  drift  velocity  of  the  projectiles  as  well  as  ionization  and  excitation  rates  of 
the  background  gas,  etc.).  The  equilibrium  values  of  the  swarm  parameters  depend  only  on  the 
ratio  E/N,  the  ratio  of  the  applied  external  electric  field  to  the  gas  number  density.  In  order  to 
assess  the  accuracy  of  the  present  numerical  results  in  the  comparative  study,  we  first  calculated 
drift  velocity  of  electrons  in  pure  gaseous  neon  for  which  considerable  experimental  data  is 
available  for  comparison.  Fig.  9  shows  the  drift  velocity  of  electrons  in  neon  for  several  values  of 
E/N  and  it  seems  to  agree  quite  favorably  with  the  corresponding  experimental  data  of  Pack  [3], 
Robertson  [4]  and  Nielson  [5].  This  superb  agreement  provides  confidence  in  the  accuracy  of  our 
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numerical  results  for  various  swarm  parameters.  Now,  for  the  comparative  study,  Fig.  10 
displays,  as  an  example,  the  equilibrium  values  of  the  drift  velocity  of  electrons  and  of  positrons  in 
the  same  gas  (namely,  neon)  for  the  same  values  of  E/N.  As  expected,  the  drift  velocity  for  both 
projectiles  approaches  zero  as  E/N  —  0.  However,  for  higher  values  of  E/N,  the  drift  velocity 
(and,  therefore,  mobility)  of  positrons  is  consistently  larger,  by  almost  a  factor  of  2.5,  than  the 
drift  velocity  of  electrons  in  neon.  It  is  anticipated  that  investigations  of  transport  properties  of 
several  different  projectiles  in  the  same  gas  will  assist  in  optimizing  the  choice  of  the  projectile  to 
be  used  in  the  gas  for  the  particular  application  at  hand. 

3.  Dissociation  of  vibrationally  excited  H2  bv  electron  impact 

As  part  of  our  investigations  we  have  continued  to  calculate  the  cross  sections  for  the 
processes  which  are  important  for  a  realistic  modeling  of  the  hydrogen  plasma.  In  a  typical 
hydrogen  source,  vibrationally  excited  H2  molecules  are  depleted  both  by  dissociative  electron 

attachment  and  by  pure  dissociation  of  the  molecule.  Since  previous  investigations  have  clearly 
demonstrated  that  the  cross  sections  for  dissociative  electron  attachment  to  H2  are  significantly 

enhanced  if  the  molecule  H2  is  initially  vibrationally  excited,  it  is,  then,  natural  to  ask  whether  the 
cross  sections  for  pure  dissociation  of  H2  are  also  strongly  dependent  upon  the  initial  vibrational 
state  of  H2.  In  the  present  investigations  we  have  interpreted  the  process  of  pure  dissociation  of 
H2  as  an  extension  of  the  process  of  vibrational  excitation  of  H2  by  electron  impact  in  which  the 

final  vibrational  level  lies  in  the  continuum.  In  such  a  situation  the  dissociation  of  the  molecule  will 
occur,  analogous  to  the  vibrational  excitation  of  molecule,  predominantly  via  the  formation  of  an 
electron-molecule  resonance.  The  significance  of  a  resonance  formation  for  the  vibrational 
excitation  to  occur  can  be  appreciated  by  a  simple  analogy:  A  small  bee  kicking  (or  colliding  with) 
a  large  elephant  would  not  cause  much  perturbation  to  the  elephant  unless  it  stays  and  buzzes 
around  the  big  animal  for  a  time  period  much  longer  than  the  normal  transit  flying  time  of  the  bee. 
During  a  collision  it  will  be  easier  for  an  electron  to  change  the  vibrational  mode  of  a  much  more 
massive  molecule  if  it  can  form  a  resonance  with  the  molecule.  At  low  energies  electron  can  form 
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Fig.  11:  Total  cross  sections  for  the  electron  impact  dissociation  of  vibrationally  excited  H2.  Vjis 
the  initial  vibrational  quantum  number  of  the  molecule  H2. 
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two  distinct  resonances  with  H2  molecule;  these  two  resonant  states  of  are  the  X  and  B 
states.  The  X  I.  resonance  decays  only  to  the  ground  electronic  state,  X  Z  ,  of  molecular 

u 

hydrogen  while  the  B  Z  resonance  can,  over  a  certain  range  of  intemuclear  separations,  decay 

o 

^  1  m  1  I 

both  to  the  excited  b  Z^  state  and  the  ground  X  Zg  state  of  H2.  Thus,  in  our  calculations,  we 
have  included  the  following  three  channels: 

e-  +  H2(X  ^Z^,  Vj)  ^  H2(B  2s+)  ^  e'  +  H2(b  ^Z;^)  e‘  +  H(ls)  +  H(ls)  ,  (30) 

e-  +  H2(X  ^Z^,  Vj)  ^  H2(B  -*  e-  +  H2(X-^Zp  -  e'  +  H(ls)  +  H(ls)  ,  (31) 

e-  +  H2(X  1z+  y.)  ^  H2(X  ^Z^)  ^  e'  +  -*  e'  +  H(ls)  +  H(ls)  .  (32) 

Using  the  traditional  formalism  of  the  resonance  theory,  we  have  calculated  the  sum  of  the  cross 
sections  for  the  dissociation  of  H2  by  low  energy  electron  impact  via  each  of  these  three  channels 
(Eqs.  (30-32)).  The  total  dissociation  cross  sections  of  H2  in  various  initial  vibrational  levels  Vj 
are  shown  in  fig.  1 1.  For  any  initial  vibrational  levels  Vj  the  contribution  of  channel  (31)  is 

negligible  at  all  electron  energies  while  channel  (30)  contributes  significantly  at  all  energies  where 

this  channel  is  energetically  open.  Contribution  of  channel  (32),  on  the  other  hand,  is  insignificant 
for  low  values  of  v^  but  increases  dramatically  as  v^  is  increased  and  for  large  initial  vibrational 

excitation  of  the  molecule  (vj  >  8)  the  contribution  of  channel  (32)  can  become  comparable  to  the 

contribution  of  channel  (30).  The  details  of  this  work  have  been  published  as  a  letter  in  a  refereed 
journal  [6]. 

4.  Asymptotic  form  of  the  penetration  probability  of  the  quantum  harmonic 
oscillator  into  the  classically  forbidden  region 

During  the  course  of  present  investigations,  we  came  across  an  interesting  problem  that 
apparently  has  not  been  addressed  in  any  previously  published  work.  The  problem  concerns  a 
linear  harmonic  oscillator  of  which  both  the  classical  and  the  quantum  versions  have  been  studied 
quite  extensively.  Surprisingly,  however,  little  seems  to  be  known  about  the  quantum-mechanical 
probability  P(n)  of  finding  the  oscillator  outside  the  classical  turning  points  when  the  quantum 
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number  n  is  large,  notwithstanding  the  obvious  expectation,  based  on  the  correspondence 
principle,  that  this  probability  must  somehow  tend  to  zero  as  the  quantum  number  tends  to  infinity. 
A  recent  article  [7]  presented  a  conjecture,  supported  by  a  purely  numerical  investigation,  that  this 
probability  has  the  asymptotic  form 


P(n)  =  A  -  B  +...  (33) 

as  n  —  oo,  with  A  »  0.133970  and  B  »  0.01 1907.  The  main  purpose  of  our  work  was  to  derive 
an  asymptotic  formula  for  this  penetration  probability  using  an  analytic  argument  based  on  the 
behavior  of  the  harmonic  oscillator  eigenfunctions  in  the  vicinity  of  the  classical  turning  points  for 
large  quantum  numbers.  We  have  shown  that  this  analytical  result  is  identical  in  form  to  Eq.  (33). 
We  have  derived  closed-form  expressions  for  the  coefficients  A  and  B,  and  their  values  are  found 
to  be  quite  close  to  those  quoted  above.  The  derivation  has  also  yielded  the  order  of  the  first 
neglected  term.  As  a  further  demonstration  of  the  usefulness  of  the  asymptotic  forms  of  the 
harmonic  oscillator  eigenfunctions,  a  formula  for  the  quantum  probability  density  inside  the 
classical  turning  points  in  the  large-/i  limit  is  briefly  derived  and  is  seen  to  have  a  simple  physical 
interpretation. 


The  energy  eigenstate  of  the  quantum  harmonic  oscillator  corresponding  to  the  quantum 


number  n  has  energy  = 


Hoa  relative  to  the  potential  minimum,  and  the  corresponding 


normalized  eigenfunction  is  given  by 

\\fjx)  =  (2"  nt  'Ik  exp  i-x^l2\^ 

r~h 

where  the  are  the  Hermite  polynomials,  and  x^  =  -M  — .  The  classical  turning  points  occur  at 


X  =  ±  'l2y  x^ 


where  for  convenience,  both  here  and  later,  we  put  v  =  n+  1/2.  The  probability  of  finding  the 
oscillator  outside  the  classical  turning  points  is  therefore 


Final  Technical  Report  (Grant  Number  F49620-92-J-(X)27) 


Page  28 


P(n)  =2  Tdx  l\}Ajj(x)l^  =  2(2"n!  J  dx  exp  (-x^/x^)  [H^^Cx/x^)] 


(2v  X, 


or,  changing  the  variable  of  integration  from  x  to  y  =  x/x^. 


/^x„ 


P(n)  = 


_}. _ ^ 

2"n!  Vi 


|dyexp(-y2)[H„(y)]2 


(34) 


Now,  using  the  asymptotic  properties  of  the  Hermite  polynomials  we  have  shown  that  P(n)  can  be 
expressed  as 

P(n)  =  2Ij  I3  +  1 12)  +  0(v-5/^)  (35) 


with 


I,  =  JdxAi2(x)=^[r(|j 

CX3 

I2  =  ^ 

00  00 

I2  =  J”dx  x^  Ai'(x)  Ai(x)  =  -Jdx  x  Ai^(x)  =  -I2 


where  Ai(x)  is  the  standard  Airy's  function  which  is  related  to  the  modified  Bessel  function  as 


Ai(x) 


'iVT’^i 


r^.3/2 


^1/3  [3 

The  integrals  1^  and  I2  are  easily  evaluated  using  the  standard  integral  [8], 

Jdxx“-'Kj{x)  = 

The  asymptotic  form  of  the  penetration  probability  is  thus  found  to  be  approximately 
P(n)  =  A  (n  +  1/2)-^^  -  B  (n  +  1/2)'^  +  O  ((n  +  1/2)-^/^) 

with 

A  =  ^  0.133975,  B  =  - ^  =  0.012252 

271^  L  1,3 jj  157iV3 


(36) 
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in  good  agreement  with  the  values  quoted  above  in  the  introductory  paragraph,  which  were 
obtained  by  purely  numerical  methods  [7]. 

It  is  of  interest  to  examine  the  absolute  error  in  the  first  two  terms  of  Eq.  (33).  Denote  this 
quantity  as  5  (n),  i.e.,  define 

5  (n)  =  P(n)  -  A  (ii  +  1/2)'^^  +  B  (n  +  1/2)-^ 

In  fig.  12  the  absolute  error  6(n)  is  plotted  against  (n  +  1/2)’^^^  for  n  =  50,  60,  70, ...,  350,  400, 
500  using  the  coefficients  A  and  B  given  by  Eq.  (36).  The  values  of  P(n)  were  obtained  by 
numerical  integration  of  the  right-hand  side  of  Eq.  (34)  with  respect  to  y^,  using  Simpson's  rule, 
on  a  grid  of  step  size  0.01,  in  REAL*8  precision  on  a  SGI  Iris  Indigo  computer.  The  solid  line  in 
the  figure  is  the  best  linear  fit  through  all  the  points.  It  is  noteworthy  that  all  of  the  values  of  5(n)  lie 

very  nearly  on  this  line,  a  fact  which  is  consistent  with  the  0-term  in  Eq.  (35). 

It  is  also  interesting  to  consider  the  probability  density  in  the  state  inside  the  classical 

turning  points.  In  elementary  textbooks  on  quantum  mechanics  [9, 10]  one  often  sees  plots  like 
the  one  in  fig.  13,  showing,  for  some  large  value  of  n  (in  this  case  n  =  50),  the  probability  density 
as  calculated  from  quantum  mechanics,  compared  with  the  corresponding  classical  probability 
density 

Pgj(x)  =  (2  V  Xq  -  x^)'^^ 

for  a  harmonic  oscillator  of  the  same  energy.  The  classical  probability  density  is  easily  derived  by 

considering  the  fraction  of  time  the  oscillator  spends  in  the  interval  between  x  and  x  +  dx;  that 
fraction  is  just  Pgj(x)  dx.  Intuitively  one  sees  that  the  "envelope"  of  the  quantum  curve  follows  the 

classical  curve  more  and  more  closely  as  n  increases.  This  can  be  seen  analytically  as  well;  yet  this 
fact  is  not  shown  in  any  textbook  or  paper  of  which  we  are  aware.  Using  an  asymptotic  series 
which  describes  the  behavior  of  the  Hermite  polynomials  for  large  order  we  have  shown  that 
l¥n(x)l^  =  P^j(x)  |2  sin^  ’  I  ®  ’ 
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Fig.  12:  The  absolute  error  6(n)  for  n  =  50, 60, 70, ...,  350, 400,  500,  and  the  best  straight  line  fit 
through  all  the  points. 
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Fig.  13:  The  dimensionless  probability  density  x^l\i/j^(x)l^  of  a  quantum  oscillator  for  n  =50  and  the 
corresponding  probability  density  (2v  -  for  a  classical  oscillator  of  the  same  energy. 
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as  n  —  oo.  Thus  to  leading  order,  the  quantum  probability  density  inside  the  classical  turning 
points  is  indeed,  for  large  quantum  numbers,  a  rapidly  oscillating  function  modulated  by  a  more 
slowly  varying  "envelope"  consisting  of  its  classical  counterpart. 

In  conclusion,  we  have  analytically  derived  an  asymptotic  formula  for  the  penetration 
probability  of  the  quantum  harmonic  oscillator  into  the  classically  forbidden  region  for  large 
quantum  numbers,  and  seen  that  it  is  in  reasonably  good  agreement  with  an  asymptotic  formula 
obtained  previously  [7]  by  numerical  methods.  We  have  also  used  a  numerical  calculation  to 
support  the  derived  order  of  the  error  in  the  first  two  terms  of  this  formula.  Finally,  we  have 
shown  that  the  leading-order  term  in  the  asymptotic  expansion,  for  large  quantum  numbers,  of  the 
quantum  probability  density  inside  the  classical  turning  points  has  a  form  which  is  readily 
understood  physically  in  light  of  the  correspondence  principle.  The  details  of  the  analytical  work 
are  in  our  paper  which  has  appeared  in  the  American  Journal  of  Physics  [11]. 

5.  Presentations  /  Publications 

During  the  tenure  of  the  present  Grant  we  have  made  presentations  of  our  research 
investigations  at  various  national  and  international  scientific  conferences.  These  are 

•  "Time-dependence  of  swarm  parameters  of  electrons  in  argon",  (with  A.  A.  Sebastian); 

presented  at  the  44th  Annual  Gaseous  Electronics  Conference,  Albuquerque,  New  Mexico, 
October  22-25, 1991. 

•  "The  mass  dependence  of  cross  sections  for  vibrational  excitation  of  diatomic  molecules  by 

electron  impact",  (with  D.  E.  Atems);  presented  at  the  44th  Annual  Gaseous  Electronics 
Conference,  Albuquerque,  New  Mexico,  October  22-25, 1991. 

•  "A  new  isotope  scaling  law  for  vibrational  excitation  of  diatomic  molecules  by  electron 

impact".  Bull.  Am  Phys.  Soc.  32,  1093  (1992)  (with  D.  E.  Atems);  presented  at  the  1992 
Annual  Meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics,  Chicago, 
Illinois,  May  19-22, 1992. 
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•  "Electron  transport  in  argon  and  neon  under  the  influence  of  an  rf  electric  field",  Bull.  Am 

Phys.  Soc.  22,  1114  (1992)  (with  A.  A.  Sebastian  and  P.  J.  Drallos);  presented  at  the  1992 
Annual  Meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics,  Chicago, 
Illinois,  May  19-22,  1992. 

•  "Transport  of  electrons  in  hydrocarbon  gases  at  high  E/N",  (with  A.  A.  Sebastian  and  P.  J. 

Drallos);  presented  at  the  45th  Annual  Gaseous  Electronics  Conference,  Boston, 
Massachusetts,  October  27-30, 1992. 

•  "Effect  of  an  rf  electric  field  on  the  charged  particle  transport  in  rare  gases".  Bull.  Am.  Phys. 

Soc.  22,  1544  (1992)  (with  A.  A.  Sebastian);  presented  at  the  Thirty-Fourth  Annual 
Meeting  of  the  Division  of  Plasma  Physics,  Seattle,  Washington,  November  16-20, 1992. 

•  "Transport  of  electrons  in  simple  hydrocarbon  gases  subjected  to  external  electric  fields", 

(with  A.  A.  Sebastian);  presented  at  the  Eighteenth  International  Conference  on  the  Physics 

of  Electronic  and  Atontic  Collisions,  Aarhus,  Denmark,  July  21-27, 1993. 

•  "Effect  of  initial  vibrational  excitation  on  the  dissociation  of  H2  by  low-energy  electron 

impact",  (with  D.  E.  Atems);  presented  at  the  Eighteenth  International  Conference  on  the 
Physics  of  Electronic  and  Atomic  Collisions,  Aarhus,  Denmark,  July  21-27,  1993. 

•  "Electron  transport  in  hydrocarbon-rare  gas  mixtures",  (with  A.  A.  Sebastian);  presented  at 

the  46th  Annual  Gaseous  Electronics  Conference,  Montreal,  Quebec,  October  19-22, 1993. 

•  "Charged  particle  transport  in  mixtures  of  hydrocarbon  and  rare  gases  subjected  to  external 

fields".  Bull.  Am.  Phys.  Soc.  2S,  1994  (1993)  (with  A.  A.  Sebastian);  presented  at  the 
Thirty-Fifth  Annual  Meeting  of  the  Division  of  Plasma  Physics,  St.  Louis,  Missouri, 
November  1-5,  1993. 

•  "Electron  transport  in  mixtures  of  hydrocarbon  and  rare  gases  subjected  to  a  time-dependent 

rf  electric  field",  (with  A.  A.  Sebastian);  presented  at  the  1994  Annual  Meeting  of  the 
Division  of  Atomic,  Molecular  and  Optical  Physics,  Washington,  DC,  April  18-21,  1994. 
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During  the  tenure  of  the  present  Grant,  our  research  investigations  have  led  to  the  following 
refereed  publications: 

•  "Vibrational  excitation  of  H2  and  HCl  by  low-energy  electron  impact:  An  isotope  scaling 

law",  D.  E.  Atems  and  J.  M.  Wadehra,  Chem.  Phys.  Letts.  197.  525  (1992). 

•  "A  new  representation  of  the  generalized  oscillator  strength  for  bound  -  free  transitions  in 
hydrogenlike  systems",  /.  M.  Wadehra  and  S.  P.  Khare,  Phys.  Letts.  A172.  433  (1993). 

•  "Resonant  contributions  to  dissociation  of  H2  by  low-energy  electron  impact",  D.  E.  Atems 
and  J.  M.  Wadehra,  J.  Phys.  B  Letter  2^,  L759  (1993). 

•  "Asymptotic  form  of  the  penetration  probability  of  the  quantum  harmonic  oscillator  into  the 
classically  forbidden  region",  D.  E.  Atems  and  J.  M.  Wadehra,  Am.  J.  Phys.  443 
(1995). 

•  "Electron  attachment  to  molecules  at  low  electron  energies",  A.  Chu^ian,  A.  Garscadden 
and  J.  M.  Wadehra,  Phys.  Rep.  (accepted,  1995). 

In  addition,  the  following  papers  are  under  preparation  and  will  be  submitted  to  refereed  journals 
for  publication  soon: 

•  "Exact  numerical  solutions  of  the  Boltzmann  equation  for  electrons  in  gases  subjected  to 
external  electric  fields",  A.  A.  Sebastian  and  J.  M.  Wadehra,  (1996). 

•  "Exactly  solvable  model  for  the  time-dependence  of  charged  particle  swarms  in  gases 
subjected  to  external  fields",  A.  A.  Sebastian  and  J.  M.  Wadehra,  (1996). 
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